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I. INTRODUCTION 



The smallness of the neutrino mass is naturally explained in the seesaw 



12j with right-handed neutrinos and in the Zee model with radiative 



Recent neutrino oscillation experiments reveal that neutrinos have small but finite mass 
differences ODD . The precise observation of the cosmic microwave background (CMB) 
shows that not only the mass differences, but also the masses themselves, are small % 

mechanism |l£J, I 1 11. 1 

mass generation |l3|. In the present paper, we focus on the Zee model, which has many 
predictions for low energy physics. This model has two Higgs doublets and a singlet Higgs 
scalar field (the Zee singlet), in addition to the fields in the standard model. In the original 
Zee model, both of the two Higgs doublets couple to the ordinary matter fields. Wolfenstein, 
however, imposed a discrete symmetry on the Zee model so that only one of the two Higgs 
doublets couples to the ordinary matter fields. As a result, dangerous processes through 
the flavor changing neutral current (FCNCs ) are easily avoided in this model, and this 
simplified model has been mainly considered so far Q|. We call this model the restricted 



Zee model (RZM) |15| . One of the attractive points of the Zee model in this version is 
that the bilarge mixing angle pattern, with respect to the solar and atmospheric neutrino 
oscillations, is naturally explained. This attractive feature, ironically, has forced the Zee 
model into a difficult position very recently. The predicted large angles must be just around 



and are incompatible with the large mixing angle (LMA) Mikheyev-Smirnov- 



Wolfenstein (MSW) solution with 6 Q ~ it/ 5.6, while the KamLAND experiment reports that 



only the LMA MSW solution is allowed [3j]. Severa 



this difficulty: introducing a fourth neutrino 



19 



ideas have been proposed to escape from 



scalar field 



2l|, introducing triplet Higgs scalars (22 



2C ] . introducing a doubly charged Higgs 



23 



24j . and allowing both of the two 



Higgs doublets to couple to the ordinary matter fields. 

Aside from the difficulty concerning the solar mixing angle, we have another issue with this 
model in cosmology. Since the lepton number is violated in this model, the primordial barvon 



number is washed out through the sphaleron process at quite a low temperature 
This model therefore seems to be unable to exist together with baryogenesis scenarios in 
which the baryon number is produced at a high temperature. However, it has been pointed 
out that there exists an approximately conserved number L' = L e — — L T for the low 
mass, low probability (LOW) and vacuum oscillation (VO) solutions of the solar neutrino 
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problem in the RZM, and complete baryon number washing out does not occur (27 1. 

In the present paper, we focus on the original Zee model in which both of the two 
Higgs doublets couple to the ordinary matter fields. We call this model the general Zee 
model (GZM) in order to clearly distinguish it from the RZM [l^ ]. Although this model was 
proposed earlier, it has not been extensively discussed, as we mentioned above. We shed light 
on the GZM once again and find not only that the two phenomenological problems above 
can be avoided j3] but also that the GZM can be consistent with experiments involving 
FCNCs. This model is attractive in the sense that the extension from the RZM, of which 
properties have been extensively discussed, is minimal. Although the extension is minimal, 
many undetermined parameters arise and the investigation in the full parameter region is 
very difficult. We then constrain the structure of the neutrino mass matrix obtained from 
the Zee model. Two constraints for the neutrino mass matrix have been proposed: imposing 
a texture to have two independent zeros in the base, where the masses of the charged leptons 
are diagonalized j^jj, and requiring the mass matrix to have a vanishing determinant 3^ |. 
While the vanishing determinant constraint is good in that the condition does not depend 
on the weak base we take, and it is suitable for an investigation into the Zee model with a 
doubly chafed H lggS si„ g let HQ, it doe 8 not fit the structure of the maS8 matrix obtained 
from the Zee model considered here. We therefore employ the two-independent-zero-texture 
constraint, of which phenomenological aspects were investigated in Refs. 

The 

suitable two-zero-texture mass matrix for the Zee model considered here is the case C defined 
in Ref. H, 

X X X ^ 

M„ = x x . (1) 
^ x x j 

The main aim of the present paper is to search for the simplest GZM with this structure 
of the mass matrix, which is compatible with the LMA MSW solution of the solar neutrino 
problem and free from the baryon number washing out J^. It turns out that such a GZM 
can be constructed by minimal extension from the RZM, i.e., by adding only one coupling 
constant to the RZM. We also examine whether this model is allowed from the viewpoint 
of experiments: /i(t) decay rate, ji{r) — > jl^] and neutrinoless double beta decay j^. 
We briefly consider the constraints from the recent CMB observation jj, 0, (], Q| and the 
possibility of the Z-burst scenario in this model [di! 



The outline of this paper is as follows: In Sec. II, we first analyze the structure of 
the neutrino mass matrix that is compatible with SNO/S-K and CHOOZ data. We next 
derive the Majorana neutrino mass matrix that is induced through the radiative correction 
in the GZM. The parameters of the GZM are determined by comparing the structure of the 
neutrino mass matrix obtained from the experiments with that obtained from the GZM. 
In Sec. Ill, we investigate the cosmological implication of the GZM. From the viewpoint 
of baryogenesis, there must be some combinations of the lepton numbers that avoid the 
washing out of the primordial baryon number. We first examine what combinations of 
lepton numbers are conserved in each type of GZM. We next obtain a range of parameters 
in which the baryon number is not washed out. We finally calculate the resulting lepton 
and baryon number. Section IV is devoted to a summary. 



II. THE GZM MASS MATRIX COMPATIBLE WITH THE LMA-MSW SOLU- 
TION 

In the present section, we review the structure of the neutrino mass matrix arising in the 
GZM and determine the parameters of the GZM using the neutrino oscillation data. We 
further examine whether the GZM is consistent with the constraints from other phenomeno- 
logical experiments. 



A. Neutrino mass matrix that satisfies SNO/S-K and CHOOZ experimental data 



The Majorana neutrino mass term is written as 



£ Majorana = -V L M U V L + H.C. for V L 



(2) 



where M. v is a complex valued symmetric matrix. Here, the weak eigenstate vl is taken 
to diagonalize the mass matrix of the charged leptons. Using the Maki-Nakag awa- S akat a 
(MNS) matrix U, we can diagonalize M. v to Ai u , 



M v = U T M V U, 



(3) 
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^ = U^ L = u 2 , (4) 

where is in the neutrino mass eigenstate. We write a diagonal matrix with a caret 
hereafter. The Lagrangian density is written in this base, 



^Majorana = ~ L ) C M u l^ + H.C. 



(5) 



We assume that M. v is real and the MNS matrix is orthogonal for simplicity. Under these 
assumptions, we can parametrize the orthogonal matrix U as 

( \ o ^ ( n «,„ \ / ft„ «-,„ n\ 

C23 S23 
\ -s 23 c 23 J 



U 



C13 S13 

1 

V -S13 ci 3 y 



C12 S12 
-S12 C12 

1, 



13l< 7 13j-<M2 



#12 (— #12), 



(6) 



where c 7 - 



cos 0jj and = sin . 



We summarize the current results of neutrino experiments to make clear what type of 
neutrino mass matrix is allowed. The S-K experiment shows that there exist a mass squared 
difference A a and a mixing angle 8 atm in order to explain the atmospheric neutrino oscillation 
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m of 



1.6 x 1(T 3 eV 2 < A a < 4.0 x 10" 3 eV 2 , 0.88 < sin 2 29 atm < 1.0 (90% C.L. 



(7) 



with the best fit values A a = 2.5 x 10~ 3 eV 2 and sin 2 29 atm = 1.00. The global analysis of 
the first results of KamLAND combined with existing data from solar neutrino experiments 
shows that there exist a mass squared difference A s and a mixing angle 9 & for explaining 



the solar neutrino oscillation 



s son 

DD 
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of 



5.1 x 10~ 5 eV 2 <A S < 9.7 x 10~ 5 eV 2 , 1.2 x 10~ 4 eV 2 < A, < 1.9 x 10~ 4 eV 2 , (8) 



0.29 < tan 2 9 Q < 0.86 (3a level), 



(9) 



with the best fit values, A s = 6.9 x 10~ 5 eV 2 and tan 2 9 Q = 0.46 |46j. This analysis concludes 
with three remarks: first, the original LMA region is narrowed and is separated into two 
islands (jSJ); second, the LMA is the unique solution to 4<r level; thirdly, the maximal solar 
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neutrino mixing angle is excluded at the 3a level Q- These data are rewritten using the 
mass squared difference ratio R and the deviation of the solar neutrino angle from 7r/4, 
= 7r/4 - e 3 , as 



1.3 x 1(T 2 < R = < 1.2 X 10 -1 , 0.038 < e 3 < 0.29 (3a level), 



(10) 



with the best fit values R = 2.8 x 10 2 and e 3 = 0.19. The CHOOZ experiment shows the 
upper bounds for the mixing angle 0i 3 (^tJ to be 



sin 2 201, < 0.1 



(11) 



Since we assume that the mass matrix M v is real, we can generally parametrize it as 

/ „ ±.h _h\ / o i i \ /oo ^ 



a +b —b 
-bed 



+ 7] 



^—b d c j 
= Ml + Ml + Ml 



1 
VI °/ 



1 

V° -iy 



Using the orthogonal matrix 



U = R2 3 {--)Ri3(0)Ri2{-^ + e 3 ) with tan2e 3 



c — a — d 



the first matrix M® in Eq. ()13|) is diagonalized as follows Q|: 

i -V2b + \{a + c - d) \{a-c + d) 



2y/2b 



(12) 
(13) 

(14) 



U£M u U = Ri 2 (-e 3 

( a+c~d+^/ (a-c+d) 2 +8b 2 



\ 



l(a — c + d) 




V2b + 



d) 







i*i2(e 3 ) 



2 









a+c-d-y/(a-c+d) 2 +8b 2 
2 



V 



(15) 



c + d J 

\ 




c + d 

According to the data ((Zj) and (jlljl . the deviations of the angles 023 and 0i 3 from the values 
023 = vr/4 and 0i 3 = are small, and the influence of r] and £ is expected to be treated 
as a small perturbation. Defining these deviations as ei = ir/A — 823 and 62 = 0i 3 , we find 
that they should lie in the range < €\ < 0.18 and < e 2 < 0.16. We obtain the following 
relations within the first order approximation: 



a — (d + c) 
£ = -2dei - V2be 2 . 



(16) 
(17) 
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Judging from the range of e\ and €2, we find that 77 and £ are much smaller than the elements 
of .M°. We therefore assume that 77 and £ are zero in the present paper. We further assume 
c = for simplicity. In this case, the structure of the mass matrix belongs to the type-C 
texture with two independent zeros, defined in Ref. j^], and the phenomenology is discussed 



in Refs. 135 



36 



22|- 



We obtain the three mass eigenvalues and e 3 as follows: 



/ a-d+^/{a+d) 2 +%b 2 



M, 



~ U^M u Uo 



2 








a-d- (a+d) 2 +8b 2 



V 



2 





( mi ^ 

m 2 
V m 3 J 



tan2e? 



a + d 



2y/2b' 

Since we identify A a and A s as A a 
the ratio R, 



18) 



(19) 



\m\ — mgl and A s 



R 



As 
A a 



4 


(a-d)^J(a + d) 2 + 8b 2 




[a - d + ^(a + d) 2 + 8b 2 


2 -Ad 2 



I m 2 — m 2 1 , we obtain 



(20) 



It is shown in Ref. |2l( that a must have an appropriate magnitude obtained from the 
experimental data mentioned above. This can be understood as follows. If we set a = 0, we 
obtain a simple relation between tan 2e% and R, 



tan 2e,q 



2V2 



4 



-R, 



(21) 



where b ~ — JA a /2 and d ~ Iv^A^-R. On the other hand, R and tan2e 3 must satisfy 



3.2 x l(T d < -R < 3.0 x 10" 2 , 7.7 x 10~ 2 < tan2e 3 < 0.66 
4 ' 3 



(22) 



from the data (jlOj) . We therefore find that a must have an appropriate magnitude because 
the conditions (J21J) and (J22J) are incompatible. Similarly, we can show that d must also have 
an appropriate magnitude. From Eqs. (|19|). (j2*U|) . and ()10j). we find that mild fine-tuning 
I a + d\ ^> I a — d\ is needed to fit the model with the LMA MSW solution. 

Because the difference a — d (= e) is small in comparison with b, a, d, it can be treated 
as a small perturbation. We obtain the following values for b, a, d, and e: 

2 

(tan 2e 3 ) 



b ~ — 1 



-3.5 x 10~ 2 eV, 

2.0 x 10~ 2 eV, 



'A a - 



(23) 
(24) 



d ~ (tan2e 3 )VA a - e ~ 1.9 x 10" 2 eV, 

= ~ 6.5 x 10~ 4 eV, 



= ~ -\ A a R 

.11 + (tan2e 3 ) 2 



for the best fit values A a = 2.5 x 1(T 3 eV 2 , e 3 = 0.19, and R = 2.8 x 10" 



(25) 
(26) 



B. GZM neutrino mass matrix 

In the present section, we review the GZM HQ 
following: 



. The interactions of the GZM are the 



£ Yukawa 
SM 



pYukawa 
*-Zee 

p cubic 

Higgs 



lL,a^i^i, a e R + H.c, 



i=l,2 



hJMabflt^h +H.c, 

\<S>l a i((T2)ab®2,bh~ + H.C. 



(27) 

(28) 
(29) 



Here, h~ is the Zee singlet, II is the lepton doublet in a weak eigenstate, $i and $2 are the 
two Higgs doublets, T\ and T2 are two Yukawa coupling matrices, / is an antisymmetric 
coupling matrix, A is a cubic coupling constant, and a, b are indices of the SU(2)l doublets. 

We write the vacuum expectation value (VEV) of each Higgs doublet as (0°) = Vi/y/2 (i 
= 1,2). Rotating $1 and $2 as follows: 

(30) 




cos (3 — sin (3 
sin (3 cos (3 




where tan/5 = V\jv^ we make only one of the two Higgs doublets have the VEV, i.e., 
= and ((fi' 2 ) = v/y/2 = \J(v\ + 1' 2 )/2. While 4>^ is the would-be Nambu-Goldstone 
boson, the diagonalization of the mass matrix for the remaining physical charged Higgs fields 

h + and <j><y~ goes as 

c+ cH 

°1 5 J 2 



c 



Higgs 



D h+ 



Ml 



-f 2 ( M i 

using the orthogonal transformation 



An 

V2 
'+\2 



\ 



M\ 



\ 



Ml 



s; 



(31) 




cos a — sin a 
sin a cos a 




(32) 
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where the rotational angle a and the two mass eigenvalues Mi, M 2 are determined as 

-V2\v 



tan 2a 



(M lj2 



(m; + ) 2 -m^ 



M 2 + (M; + ) 2 ± J[Af* - (M; + ) 2 ] 2 + 2AV 



(33) 
(34) 



Using $1 and $ 2 m Eq. (J5U|) . the Lagrangian density (|27jl is written as 



p Yukawa 
*-SM 



Since only 02° nas the VEV, using the biunitary transformation 



<<L.a 



cos/3 r\ - sin/3 r 2 )$i Q + (sin/3 r x + cos/3 T 2 )$ 2a e fl + H.c. (35) 



( v. \ 



^L.a = V Lh,a, 1™1 = 



lm — Ttn 
) L L,2 — e L 



(36) 



The matrix (sin/3 Ti + cos/3 r 2 ) is diagonalized as 



(37) 



f m ee \/ L (sin/3 r x + cos/3 r 2 )Vj ee — M l 



l—m 



V2 

v 



m P 



m M 



(38) 



y m T J 

where Mi- m is the mass matrix of the charged leptons. The Lagrangian density (|27j) is 
written as 



nYukawa 7*" 1 

*-SM — l L,a 



(39) 



^cos/3rr-sin/3r-)$ lia + r m $ 2 j 

where 17 ee V^V^ (i = 1,2). 

Since we estimate the interaction rates in the early Universe in Sec. IIII[ we write down 
the Lagrangian density in the symmetric phase (v = 0) in addition to that in the symmetry 
broken phase (v = 246GeV) where the neutrinos have the Majorana masses. 

1. Symmetric phase (v = 0) 

The interaction Lagrangian density is as follows: 



("Yukawa / 

*-SM — L L,a 



-dm r 



"tan/3 sin/3 

{Yukawa = ^jm^ny^ + H.C, 
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(40) 
(41) 
(42) 



Here, all leptons are massless and f m (= VifV^) is antisymmetric. 



2. Symmetry broken phase (v = 246 GeV) 

The interaction Lagrangian density is as follows: 



n Yukawa 
*-SM 



fYukawa 
'-Zee 



'^-f m sin aSt + cos aS+) + f "> 2 + 

Tan/3 sm/3 2 yv 1 11 2 



e R 



-*m 



1 



,17 + (T + iiV2, 



tan/3* sin/5* z /iri v y/2 

2vZf m (?£) c (cosaS± + sinaSj) + H.c, 

r 1 



eS + H.c, (43) 
(44) 



£™f g g S = A(cosaS' 1 +sinaS' 2 

H.c. 



y/2 



[a + iN 2 ) (— sin aSf + cos aS^ 



>;v; 4 



(45) 



Here, the Higgs fields S 1 , 5 2 have masses Mi, M 2 , respectively. The fields </> 2 + an d -^2 
are the would-be Nambu-Goldstone bosons. 

In this phase, the Majorana neutrino mass is calculated from the two one-loop diagrams 
in Fig. ^and the transposed diagrams of them, for m e , m M , m T -C Mi, M 2 , 

a/3 



iM 



a 6 



A 



{fMf + Mff 



^/2cos/5 



(46) 



iM 



where A = (l/8V2ir 2 v tan /3) sin 2a log(Mf/Mf) and 

72 ^ ) 7°^ are elements of the matrices 

r™, r m , respectively. 



Si 



— * * <— 



iM 



el 



J_Y*__Ly 2 * 

yt R s a 



Py 



a 

Vl 



el 



'lis f<°< 



— ic 



t 

1 - 



-J-y: 



FIG. 1: The diagrams that are used to calculate the radiatively induced neutrino mass matrix. We 
use s a = sina,c a = cosq, and t a = tana and so on in this figure. 
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C. Matching the GZM with the experiments 



We obtained the relation between the value of a and the experimental data in Sec. Ill Al 

a = tan2e 3 A/A^. (47) 



On the other hand, we calculated the radiatively induced a in the GZM in Sec. Ill Bl 

a = -2AB[m,r^*f + m./'^T], (48) 

where B = v/y/2 cos/3. From Eq. (|37|) and Eq. (|4*%|) . we find that one of 7f T and 7^ 
must be nonzero at least. Because we prefer a minimal extension here, we assume that the 
other elements of T™ are zero. From now on, we further assume that / and r™ are real. 
Here, we match the neutrino mass matrix in the GZM (jUaj) with that determined from the 
experimental data (J2*3|) — as follows 

b = -y^ = A{rml - Bm T r^} = -A\f^m 2 T + Bm^T], (49) 
a = tan2e 3 = -2AB[mJ e ^T + m T f eT ^ T ], (50) 
d = tan2e 3v // Al-e = Af^ T m 2 T . (51) 

In this model, there are nine parameters (Mi, M2, a, (3, f efl , f eT , f^, 7f r , 72 M )- From Eqs. 
(}4T?j) — ()51[) . we parametrize four parameters (/ eM , / eT , / Mr ,72 r ) using the other five param- 
eters (Mi, M 2 , a, /3, 7^) as follows 

1 t[6 + Bj(m>g)2f1 + 
J ' AmJ[6 + Bd(m„/m?)l5'']-.Bdm,/y?'' 1 ' 

- ^ < 54 > 

= m T b^f + ( y a/2B)m fl 
72 mpb + Bd-yFil/m* -1/m 2 )' 

We consider the simplest two cases, (1)72 T 7^ 0,7^ = and (2) 7| T = 0, 7^ 7^ 0. In these 
cases, we have the following relations. 

(1) 'y| T ^ 0,7^ = 0. Equations (O - (gSJ) are written as 

; " Aml + 2Abml^ 6< * iU ^ A J' 

11 



d 8 /lO- 4 GeV _1 \ PT am r ,/ cos 3 \ 

f = ~ 6.1 x 1(T 8 - , 7f = — f ~ -2.1 x 1(T 3 . 

J Am 2 T \ A y 12 2Bb \VV2j 

(2) 7| r = 0, 7^ 7^ 0. Equations (J52j) — are written as 

F T = -A * 6.1 x 10~ 8 ( 10 " 4 G A eV " ) , 7f = - 1.2 x 10- ( ^l) . 

J Ami \ A J 2Bb \VV2J 

D. Experimental constraints on r(/x) decay, r{p) — > e 7 and Oz>/3/3 amplitude 

In the present section, we examine whether the parameters obtained in the previous 
section [Eqs. ()5bj) .(|57| ) ] are consistent with the constraints from the processes r(/i) decay, 
r(yu) — > e7, and 0i/38 decay jl^ ]. 

1. r(/i) decay 

Since the leptonic flavor changing neutral current does not exist in the standard model 
at all, stringent constraints are imposed on their decay amplitudes 3]. There is, how- 
ever, r(p) decay through the FCNCs at the tree level in the GZM. These constraints 
are satisfied for the parameter set (jHfij) and (JH7jl as follows for 3 ~ 71/ 4: 

1 ^ ' G^'i) 4 n 

« 1.8 x 1(T 6 (experiment |48(), 

Br(/i" e-e-/i + ) ~ ' 2~,„1J ~ 8 - 7 x 10 



<C 1.0 x 10 12 (experiment 



48|), (58) 



where Gp is the Fermi coupling constant and M'\ ^100 GeV. Since the coupling 
constants / eM , f eT , and f^ T are much smaller tha n 7o T a nd 7^, the contributions to 
the branching ratios purely from / are negligible j^iil Isoi Isij l . 

2. r(/i) -> e 7 
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Since these processes violate the lepton flavor, they do not exist in the standard model 
at all. On the other hand, these processes become possible in the GZM through the 
exchange of the charged and neutral Higgs bosons. Since they are rare events, stringent 
constraints are imposed on the amplitudes j^. These constraints are safely satisfied 
for the parameter set (|56|) and (|57|) as follows for (3 ~ it /A: 



\Y 2 T (V2m T /v)\ 



2 



Bi(r- e"7) ~ a^-^ < 2.3 x KT 10 < 2.7 x 10~ 6 (experiment |4 



GUM'i) 
h efl (\/2m /v)\ 2 

Br(u- -> e~j) ~ a^^Ar tt^ £ 2.8 x 10^ 15 < 1.2 x 1CT 11 (experiment |48jl). 

; G^(M'?) 4 v ^ 

(59) 

Again, since the coupling constants / eM , / eT , and f^ T are much smaller than t^T and 
72 M , the contributions to the branching ratios purely from / are neglig ible QMS- 

3. Oupp 

The observation of the neutrinoless double beta decay (OvfiP) gives evidence that 
the neutrinos have Majorana masses. The study of this process is a crucial test of the 
validity of our model, since neutrinos have Majorana masses, the interaction of the Zee 
singlet h~ violates lepton number, and the Yukawa coupling is flavor changing. As 
was discussed by Schechter and Valle ^ in the theory with an exotic doubly charged 
scalar, the extension of the Higgs sector possibly provides a new source of the decay. 
We can, in fact, think of a process with the exchange of our exotic scalar, i.e., h~ [see 
Fig. EK a )]- However, we readily know that this diagram does not exist in our model. 
This is because the Yukawa coupling of h~ is flavor off-diagonal and the intermediate 
neutrinos are either v ^ or v T . This, in turn, means that the Yukawa coupling of 
<pi at another vertex should be flavor changing, i.e., 7^ or 7™, which vanishes in our 
model. Another source of the double beta decay comes from the exchange of Majorana 
neutrinos with Majorana mass insertion due to W ± exchange. As was first pointed 
out by Wolfenstein j^j, the amplitude is proportional to m ee . The absence of any 
report of this process so far gives the upper bounds to the element m ee in the neutrino 



mass matrix 



60], 



< 0.1 eV (90% C.L.). (60) 
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The value obtained in Sec. Ill Al is consistent with this upper bound as follows: 

m ee = a ~ 2.0 x 10~ 2 eV. (61) 

Since the Yukawa coupling is flavor changing, another source of decay exists as 
shown in Fig. |2fb), where m wt or m TT contributes, but not m ee . This diagram, 
however, does not exist from reasoning similar to the above; namely, the amplitude of 
the diagram is also proportional to or 72 6 , which does not exist. Thus the ordinary 
process due to m ee discussed above is the unique source of double beta decay. If the 
sensitivity of the experiment is improved by a factor of 5 or so, the decay should be 
seen if our model is valid. 




v 



(b) 



FIG. 2: The diagrams contributing to neutrinoless double beta decay. In (a), the h~ interaction 
violates the lepton number, while in (b) the lepton number is violated by the insertion of m^. 
Both processes fail to exist because of the vanishing Yukawa coupling 72 6 = 0. In both diagrams 
Uft may be replaced by u T . 



III. COSMOLOGICAL IMPLICATIONS OF THE GZM 

In Sec. HH we obtained the values for the coupling constants (/ eM , f eT , /^ r , 7| T , 7^) sug- 
gested by the neutrino oscillation data. In the present section, we investigate whether the 
primordial baryon number is washed out or not for these values. If there are some leptonic 
conserved numbers, the primordial baryon number remains finite even though the sphaleron 
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process partially washes out the baryon number. We therefore consider what leptonic num- 
bers are effectively conserved at a high temperature. We then calculate what amount of the 
primordial baryon number remains in the case that there exist some conserved numbers. We 
also briefly consider the constraints from the recent CMB observation 

110,0,0 and the 



possibility of the Z-burst scenario in this model |4ol. |41|. 



A. Conserved number for each coupling 

It is convenient to define the following global U(l) charges as 

Lefj, = L T — L e — L M , L eT = — L e — L T , = L e — — L T (62) 

for an investigation into the washing out of the baryon number. These quantities are con- 
served numbers in the standard model. We, however, have coupling constants which violate 
these conserved numbers in our model. Each coupling conserves two of the three indepen- 
dent numbers which are shown in the Table |l] using a linear combination of L efl ,L eT , and 



Coupling constant 




r 




if 


72 


Conserved number 













TABLE I: Conserved number for each coupling constant. 



B. Interaction rates induced by each coupling 

We here estimate the interaction rates induced by the coupling constants / and 72. 



I-IT 



1. L^ T violating process induced by the coupling constant / 

We consider the h~ decay process and the inverse decay process induced by the coupling 
constant/ MT . The interaction rates for these processes are estimated at a temperature T 

^M h yy, 

where Mh is the mass of the Zee singlet and (E) ~ 2.7 T (boson) is the averaged energy. 
Similar results are obtained for the processes induced by the coupling constants f e ^ and f er . 
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2. Processes induced by the coupling constant 7" 



We next estimate the interaction rates for the processes induced by the coupling constant 
7f T . We estimate the interaction rates of (j) t decay and the inverse decay. The interaction 
rate for the process <f> 1 



<-> [y e i) c T R is estimated for T > M-f 1- 



1 


7| T 


227T 


sin j3 



(64) 



where M ± + is the mass of the field <$> x . Similar results are obtained for the processes induced 
by the coupling constant 7^. 



C. Out-of-equilibrium condition to avoid the baryon number washing out 

In the present section, we consider which processes discussed in the previous section were 
out of equilibrium in the early Universe jsjj. We will consider what numbers are conserved 
in the case where these processes are out of equilibrium in the next section. The condition 
that the processes are out of equilibrium is given by 

r < H, (65) 

where H is the Hubble parameter, H = lM^T 2 /M Pl ~ 1.47 x l(T 18 (T 2 /GeV) with 
Mpi being the Planck mass. Here, g* is the total degrees of freedom of effectively massless 
particles. In the GZM, we adopt the value g* = gf M + g% + gf 2 ~ 112.75. 

1. Coupling constants f efl , f eT , and f^ T are out of equilibrium 

ftJ-T process. We examine the condition that the process h~ <-> is out of equilibrium 
at T = M h . We apply the condition (JoT)]) to the rate (|63|) . 

T(h- <-> yTu T ) < H (66) 

^^"4 < i - 5xio " l8 (c^)- < 67 » 

We obtain the condition for T = Mh, 



2.9 x 10 6 ^^ GeV 2 < 
sin 2a 



Ml 



[Ml + Ml + {Ml - Ml) cos 2a] , (68) 



16 



where we use the following relations: 



1/ 



i JT \ 





d 




Ami 



A 



sin 2a 



Ml 



8v^?r 2 ?;tan/3 
Ml + Mf + {M\ - Ml) cos 2a 



M 2 
l0g Mf 



(69) 
(70) 



At the higher temperatures, T > M^, the out-of-equilibrium condition Eq. ()67|) is more 
easily satisfied. At lower temperatures, T < M^, the number density of the Zee singlet h~ 
in the thermal bath is so small that the L^-violating processes are negligible. It is therefore 
sufficient to consider the out-of-equilibrium condition at T ~ M^. We set (3 = vr/4 for 
simplicity hereafter. The inequality 1)68)1 is written for this value of j3, 



2.9 x 10 6 
sin 4 2a 



GeV 2 < lot 



Ml 
Ml 



Ml + Ml + {Ml - Ml) cos 2a 



(71) 



From Eq. 1)71)1. we find that this condition is easily satisfied for the angle sin 2a ~ 1. We 
then fix the value of a to 7r/4. In this case, the inequality (|71)1 is simplified as 



18 < (log£)V + y 2 ], 



(72) 



where we parametrize as x = Mi/100 GeV, y = M 2 /100 GeV (x, y > l,x > y). From this, 
we find that this process is out of equilibrium for M\ > 310 GeV, when we fix the value of 
M 2 to 100 GeV [ej. 

jer p rocess The out-of-equilibrium condition for the f eT process is obtained just as the 
condition for the / Mr process. We examine the condition that the process h~ «-> e~v T is out 
of equilibrium at T = M h . We apply the condition to the rate (|53^l . 



5.4vr^ 
2.9 x 10 7 



M 2 

er|2 h < 1.47 x 10^ 18 



1 T 

GeV> < (l„ g §' 



sin 2a \ jh 2 / 

In the case a = 7r/4, the inequality (|74^1 is simplified as 

1.8 x 10 2 < (log-) 4 [x 2 + 



(Gev)' 

M 2 + M 2 + (M 2 - M 2 ) cos 2a 



(73) 
(74) 



y 



(75) 



From this, we find that this process is out of equilibrium for Mi > 500 GeV, when we fix 
the value of M 2 to 100 GeV. 
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jen p rocess xhe out-of-equilibrium condition for the f e>1 process is obtained in a similar 
manner to that in the f^ T process. The condition T{h~ <-> e~v^) < H is written as 

2 -° X 4o 1Ql7 GeV 2 < flog 4 [Ml + Ml + {Ml - Mf) cos 2a] . (76) 
sin 2a \ M| / L J 

From this, we find that this process is out of equilibrium for Mi > 1.2 x 10 6 GeV, when 

we fix the value of a to 7r/4 and the value of M2 to 100 GeV. Since we recognize that this 

mass range for Mi is quite large and somewhat unnatural, we consider that this process is 

in equilibrium at T = M^. 

2. Coupling constants 7f T and 7^ are out of equilibrium 

We consider which processes induced from the coupling constants 7| T and 7^ are out of 
equilibrium for each case considered in Sec. Ill CI 7f r 7^ 0,7^ = and 7| T = 0,7^ 7^ 0. 

7| r 7^ and 7^ = 0. We examine the condition that the process <$' x <-> (^£) c t^ is out of 
equilibrium at T = M h , 

T(07 ~ KTr R ) < H, (77) 

=► 8.4 x 10 10 GeV < (78) 

where we fix the value of a to tt/4. From this, we find that this process is out of equilibrium 
for Mi > 1.2 x 10 n GeV, when we fix the value of M 2 to 100 GeV. We consider that this 
process is in equilibrium at T = Mh for the same reason mentioned above. 

7| r = and 7^ 7^ 0. We examine the condition that the process 4> 1 (i^IY^r 1S ou t 
of equilibrium as follows: 

r(07 <-> (ut)^ R ) < H, (79) 

=► 2.8 x 10 8 GeV < ^ ^ , (80) 

where we fix the value of a to 7r/4. From this, we find that this process is out of equilibrium 
for Mi > 4.0 x 10 8 GeV, when we fix the value of M 2 to 100 GeV. We consider that this 
process is also in equilibrium at T = Mh- 

D. Final lepton number and baryon number 

In the present section, we consider the final lepton and baryon number for the two cases 
IT + 0, T f = and 7 2 er = 0, 7 f + 0. 



18 



1- 12 ^ 0,72^ = 



Since the processes induced by 7| r and f etM are in equilibrium from the previous section, 
the only possible conserved quantity is 

P = L eT . (81) 

This quantity is conserved if the process induced by f er is out of equilibrium [Eq. 1)74)1] . 
The region in which this quantity is conserved is shown in Fig. El This conserved number is 
converted to the baryon number through the sphaleron process and the final baryon (B^ q ) 
and lepton number {L eq ) are calculated as 

Bp = «L(* + u 

f 563^3 



L f = -—Bf, (82) 
f 20 f 



where (B/3 + L er )i is the primordial value generated at the temperature much higher than 
T = M/j. The region shown in Fig. |3]is therefore the allowed region for avoiding the washing 
out of the primordial baryon number in this case. 

2- 7f T = 0,7fV0 

Since the processes induced by 7^ and / eM are in equilibrium, the only possible conserved 
quantity in this case is 

P = L eT + hy. T {= 2L T ). (83) 

This quantity is conserved if the process induced from f eT and f^ T is out of equilibrium [Eq. 
(174)1 and Eq. 1)71)1 ]. The region in which this quantity is conserved is shown in Fig. El This 
conserved number is converted to the baryon number through the sphaleron process and the 
final baryon (B'p) and lepton number {Lj 1 ) are calculated as 

f 185 V 3 



L f = --B?, (84) 
f 20 f 



where (B/3 — L T )i is a primordial value generated at the temperature much higher than 
T = Mh- The region shown in FIG. El is therefore the allowed region for avoiding the 
washing out of the primordial baryon number in this case, also. 
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FIG. 3: The primordial baryon number remains finite in the allowed region. The out-of-equilibrium 
condition 1.8 x 10 2 < [log(y/x)] 4 [x 2 + y 2 ] [Eq. (|75j)] is used in this figure with the variables 
x = Mi/100 GeV, y = M2/IOO GeV. Only the region x > 1, y > 1, x>yis considered. 



E. WMAP and Z burst 



The investigation into the cosmic microwave background by the Wilkinson Microwave 
Anisotropy Probe (WMAP) recently gave the following upper bound to the sum of the 



neutrino masses 



Kl < 0.71 eV (95% C.L.). (85) 

i 

This value is consistent with that obtained in Sec. Ill Al 



Y Kl = vO - e ) 2 + 862 + d - °- 13 eV - ( 86 ) 

i 

If the ultrahigh energy cosmic rays with energy beyond the Greisen-Zatsepin-Kuzmin 
limit are explained in the Z-burst scenario, the heaviest neutrino mass should lie in the 



range 



are expiame 

0.08 eV < I m v (heaviest) | < 1.3 eV (68% C.L.). (87) 



The value obtained in Sec. Ill Al 



(2a - e) 2 + 86 2 



\m 1 \ = 2 ~5.3xlO~ 2 eV, (88) 
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is near this range, and the ultrahigh energy cosmic rays might be explained in the Z-burst 
scenario in our model. 



IV. SUMMARY 



We propose a simple GZM which is compatible with all the neutrino oscillation data for 
the LMA MSW solution of the solar neutrino problem, and is free from the baryon number 
washing out. Our extension from the RZM is minimal, i.e., we add only one coupling 
constant 7f T or 7^ in addition to those in the RZM. To avoid the baryon number washing 
out, the masses M\ and M2 must be in the region shown in Fig. |3J The final baryon and 
lepton numbers are calculated for the two cases (l)7| r 7^ 0, 7^ = and (2) = 0, iT ; 

(1) H / 0,-ff- = : Bf = ||(f + L„) , L? = -f B ? , (89) 

(2) 72 - = 0, 7 f + : Bf = ^(| - L T ), Lf = -|flj«. (90) 

In order to explain today's baryon number, it is therefore necessary that the quantity (-B/3 + 
L eT ) or (B/3 — L T ) is generated in the early Universe or the baryon number should be 
produced at lower temperatures T M^. We check that the predicted lepton flavor violating 
processes are not in conflict with the phenomenological constraints for r(//) decay, r(/i) — > 
ej, and 0vj3j3 decay. 
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